Astronomy & Astrophysics manuscript no. asas'astroph 


©ESO 2013 


March 4,2013 





A new bright eclipsing hot subdwarf binary from the ASAS and 

SuperWASP surveys 

V. Schaffenroth 1 - 2 , S. Geier 1 ' 3 , H. Drechsel 1 , U. Heber 1 , P. Wils 4 , R. H. 0stensen 5 , P. F. L. Maxted 6 , and G. di Scala 7 



co 
o 



6 

-t— > 



Dr. Remeis-Observatory & ECAP, Astronomical Institute, Friedrich- Alexander University Erlangen-Niirnberg, Stern wartstr. 7, 
96049 Bamberg, Germany 

e-mail: veronika. schaffenroth@sternwarte.uni-erlangen.de 

Institute for Astro- and Particle Physics, University of Innsbruck, Technikerstr. 25/8, 6020 Innsbruck, Austria 
European Southern Observatory, Karl-Schwarzschild-Str. 2, 85748 Garching, Germany 
Vereniging voor Sterrenkunde, Belgium 

Institute of Astronomy, KU Leuven, Celestijnenlaan 200D, B-3001 Heverlee, Belgium 

Astrophysics Group, Keele University, Staffordshire, ST5 5BG, UK 

Carnes Hill Observatory, 34 Perisher St . , Horningsea Park, NSW, Sydney Australia 



Received 14 December 2012 / Accepted 01 March 2013 



ABSTRACT 



We report the discovery of a bright (m v = 11.6 mag) eclipsing hot subdwarf binary of spectral type B with a late main sequence 
companion from the All Sky Automated Survey (ASAS 102322-3737.0). Such systems are called HW Vir stars after the prototype. 
The lightcurve shows a grazing eclipse and a strong reflection effect. An orbital period of P = 0. 13927 d, an inclination of i = 
65.86°, a mass ratio q = 0.34, a radial velocity semiamplitude K l = 81.0 km s -1 , and other parameters are derived from a combined 
spectroscopic and photometric analysis. The short period can only be explained by a common envelope origin of the system. The 
atmospheric parameters (T e g = 28 400 K, logg = 5.60) are consistent with a core helium-burning star located on the extreme 
horizontal branch. In agreement with that we derived the most likely sdB mass to be M s(iB = 0.46Mq, close to the canonical mass 
of such objects. The companion is a late M-dwarf with a mass of M corap = 0.16Mq. ASAS 102322-3737.0 is the third brightest of 
only 12 known HW Virginis systems, which makes it an ideal target for detailed spectroscopic studies and long term photometric 
monitoring to search for period variations, e.g. caused by a substellar companion. 

Key words, stars: subdwarfs, binaries: eclipsing, binaries: spectroscopic, stars: early-type, stars: fundamental parameters, stars: 
individual: ASAS 102322-3737.0 
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Hot subdwarf stars (sdB) are evolved, compact stars found in 
the disk and halo of our Galaxy. They dominate the population 
of faint blue stars. Especially in the context of Galaxy evolution 
sdBs are important because they are believed to be the dominant 
source for the "U V upturn phenomenon" w hich is observed in 
elliptical galaxies ( Bro wn et al .1 I19971 12000h . Subdwarf B stars 
are core helium burning stars on the extreme horizo ntal branch 
■ (EHB ) . They have very thin hydrogen envelopes (Saffe r et al.l 
1 1994 iHeberl Il986l) which avoid hydro gen shell burning. 
Therefore they evolv e from the EHB directly to white dwarfs 
dDorman et al.lll993l) . 

The formation of sdBs requires an extraordinarily high mass 
loss on the red-giant branch (RGB). About half of the sdB stars 
are f ound in close binari e s with periods from ~ 0.05 to 30 d 
(e.g. iMaxted et all |200il iNapiwotzki et all l2004t iGeier et al l 
120121: IVennes et al.l 1201 2h . Hence, mass transfer must play 
an important role in the formation o f thes e stars . Performing 
binary evolution studies lHan et all d2002l 120031) found that 
common-envelope evolution, resulting from dynamically unsta- 
ble mass transfer near the tip of the RGB, should produce such 
short-period binaries (P = 0.1 - 10 d). The most probable mass 
for such sdBs is 0.47 M© (lHan et al]l2003l) . which is called the 
canonical mass. 

Eclipsing post-common envelope binaries that consist of sdBs 



and late M star companions with periods of about 2 - 3 hr are 
called HW Vir systems after the prototype. Such systems are 
of high value because it is possible to derive the mass of the 
sdB as well as the mass and the nature of the companion from a 
combined analysis of time resolved spectra and the lightcurve. 
These systems are rare but can easily be recognized by the 
prominent reflection effect which is the only contribution of 
the companion to the optical light. Until now only twelve such 
systems are known. Interest in HW Virginis systems has been 
revived by the disco very of substella r companions t o the pro- 
totyp e HW Virginis dLee et al J 120091) and HS0705 (lOian et al.l 
2009 ) via the light travel time me thod (for an up-to-date census 
see IZorotovic & Schreiberl (120121) ). For this method a long term 
photometric monitoring for several years is necessary. 
ASAS 102322-3737.0 (ASAS 10232 for short) was discovered 
as a variable s t ar in course of the All Sky Automated Survey 
by iPoimanskH d2003l) . but misclassified as a 6 Scuti star. In 
2007 this system was found in the ASAS Survey by P. W. and 
recognized as a HW Virginis system. Another lightcurve was 
obtained at the Carnes Hill Observatory, Sydney, in March 
2008. BVRI photometry was taken with the Euler telescope on 
La Silla in April 2008, c onfirming the presenc e of an eclipse. In 
the SuperWASP survey dPollacco et al. 1 120061) this system was 
also observed in several runs. ASAS 10232 is the third brightest 
HW Vir system with a my = 11.6 mag, which makes it well 
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suited for follow-up observations. 



2. Observations 

2.1. Photometry 

Observations of AS AS 10232 were taken in the nights of the 
10., 11., 12., 13., and 31.3.2008 in the B, V, and I band 
with a 12" LX200 GPS Schmidt-Cassegrain Telescope at the 
Carnes Hill Observatory in Sydney, Australia. More lightcurves 
in BVRI were taken with the 1.2-m Leonhard Euler Telescope 
at the La Silla Obervatory, Chile, in the nights of the 23. and 
25.4.2008. Moreover, a lightcurve of this system was observed 
with the robotic SuperWASP telescope at the South African 
Astronomical Observatory in three runs from May 2006 to 
January 2009. For this observation a broad band filter (400-800 
nm) was used. 

2.2. Spectroscopy 

15 medium resolution spectra (R = 3400) of AS AS 10232 with 
a exposure time of 70 seconds were obtained with the EMMI 
spectrograph at the ESO-NTT Telescope in La Silla, Chile, from 
the 11. to 14.1.2008. These spectra cover only a small part of 
the optical spectrum (3900 - 4360 A) and are well suited for 
radial velocity work, but not for the determination of the at- 
mospheric parameters. Hence, another 33 spectra (R ~ 1750, 
3800 - 6540 A) were obtained on the 20. and 25.2.2010 with 
the GMOS spectrograph mounted at the 8.1 m Gemini South 
telescope at Cerro Tololo Observatory in service mode with 
an exposure time between 100 and 300 seconds. The EMMI 
data were reduced with the MIDAS package distributed by the 
European Southern Observatory (ESO). The GMOS data were 
reduced with the PAMELA^ and MOLLY 1 packages. 
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Fig. 1. Radial velocity plotted against orbital phase of 
ASAS 10232. The RV data was folded with the period measured 
from the SuperWASP lightcurve. The RVs were determined 
from ESO-NTT/EMMI spectra (triangle) and Gemini/GMOS 
spectra (rectangular). The errors are formal 1 cr uncertainties. 
The GMOS-RVs were shifted systematically to fit to the EMMI 
results. 



determined to K\ — 8 1 .0 + 3 km s _1 . Fig. Q]shows a phased RV- 
curve where the radial velocities derived from the GMOS spectra 
have been shifted by -27 km s _1 . 



3. Spectroscopic analysis 

The optical EMMI and GMOS spectra give a good phase cover- 
age of the radial velocity curve. The GMOS spectra enable us to 
determine the atmospheric parameters covering about half of the 
orbit. 



3.2. Atmospheric parameters 

Atmospheric parameters were determined by fitting synthetic 
spectra to the observed Balmer and helium lines of each of 
the 33 GMOS spectra using SPAS dHirschl 120091 A model 
grid of synthetic spectra was calculated by using LTE model 
atmospheres with solar metalicity and metal line blanketing 
dHeber et al.ll2000t) . As some HW Vir systems showed an appar- 
ent c hange of the atmospheric parameters over the orbital phase 
(e.g. lWood & Safferlll999l:lDrechsel et~alll200ll) all spectra were 
fitted separately. This effect is linked to the reflection effect as 
the contribution of the companion to the spectra is stronger the 
more the heated side is visible. 

The variation of the parameters over the phase is clearly visible 
for our star system likewise. The temperature seems to change 
about 1500 K and the surface gravity about 0.09 dex over the 
phase as it can be seen in Fig. [3] The helium abundance shows 
a scatter of about 0.5 dex, which is not much more than the sta- 
tistical error. For a grazing eclipse the smallest contamination 
by the secondary and therefore the most accurate value is ex- 
pected for the phase 0. Unfortunately, the binary was not ob- 
served at this phase. We adopted the value closest to phase zero 
and therefore a temperature of 28400+500 K, a surface gravity of 
log g = 5.60+0.05 and a helium abundance of about 1 .8+0.2. We 
considered bigger errors than the statistical errors, which can be 
seen in Fig. [3] to account for the change of the parameters over 
1 http : //www2 . Warwick. ac .uk/fac/sci/physics/research/astlte/phQ^^ftAsyblesrifttiBaaaa'ors. The helium abundance is sub- 



3. 1 . Radial velocity curve 

The radial velocities were measured by fitting a combination of 
Gaussians, Lorentzians and polynomials to the Balmer and he- 
lium lines of all spectra. Assuming a circular orbit sine curves 
were fitted to the RV data points in fine steps over a range of 
test periods. For each period the y 2 of the best fitting sine curve 
was determined (see iGeier et al Jl20TTab . At first the EMMI- 
and the GMOS-datasets were fitted separately. In each case the 
best solution had the same orbital period as derived from the 
lightcurve (P 0.139 d). Both datasets cover either the full 
phase (EMMI) or about half of the phase (GMOS). While the 
RV-semiamplitudes are similar, the system velocities of the two 
datasets show a significant shift of about 27 kms -1 with respect 
to each other. This shift is probably due to an instrumental effect. 
More observations are needed to solve this issue. 

Taking the average value of the solutions derived from the 
two datasets the semiamplitude of the radial velocity curve was 
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Table 1. Radial velocities of ASAS 10232 
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solar like in most sdBs. Such atmospheric parameters are typical 
for sdBs in HW Vir systems (see Fig.|4]i. 



4. Photometric analysis 

4.1. Ephemeris 

The SuperWASP lightcurve (see Fig. ID clearly shows that 
ASAS 10232 is a short-period binary with a grazing eclipse. The 
huge reflection effect indicates a cool main sequence star com- 
panion. The Carries Hill Observatory lightcurves (see Fig. |6]l in 
the three bands B, V and I were used to determine the time of 
the primary minimum. Parabolas were fitted to the minimum ob- 




-40 -20 20 40 

Fig. 2. Fit of the GMOS spectrum nearest to phase 0, the dashed 
line shows the measurement and the solid line shows the best 
fitting synthetic spectrum 

tained at the night of the 13.3.2008, where the minimum was 
clearly visible in all bands. From the different measured times of 
the minima in the different bands the standard deviation was cal- 
culated as the error in time. As the SuperWASP lightcurve has 
much higher accuracy and covers a longer time span it was used 
to determine the orbital period of the system. The period deriva- 
tive was found by a fit of a parabola to the O-C curve measured 
from the SuperWASP and BVI lightcurves, see Fig. [7] 

HJD = 2454538^99689 + O?13926940 ■ E - 6<?1 • 10" 11 • E 2 

+42 +4 +2.3 

The period derivative is not yet a 3 sigma detection, more 
observations are necessary to confirm it. 

4.2. Lightcurve solution 

The phased B, V, I lightcurves as well as the white light 
SuperWASP lightcur ve were analysed with the MORO code 
dDrechsel et al.ll995h . This lightcur ve solution program is base d 
on the Wilson-Devinney approach dWilson & Devinnevlll97ll) . 
but uses a modified Roche model that considers the radiative 
pressure of hot binaries. For the analysis of the SuperWASP 
lightcurve normal points were formed by binning the fluxes of 
individual measurements over narrow time intervals. This is nec- 
essary because of the huge number of data points (13816 data 
points were binned to 478 normal points) that would make the 
analysis very tedious. 

We used the Wilson-Devinney mode 2, which does not restrict 
the system configuration and links the luminosity and the tem- 
perature of the second component on the basis of the Planck 
law. As the luminosity ratio is very high and the companion con- 
tributes almost exclusively via the reflection effect, the measured 
temperature of the companion is not reliable. Because of the high 
number of parameters (12 + 5n = 32) some have to be deducted 
from spectral or theoretical constraints. 

Due to the early spectral type of the primary star the gravity 

darkening exponent can be fixed at g] = 1 as expected for 

radiative outer envelopes (Ivon Zeir)ellll924h For the cool con- 
vective companion g2 was set to 0.32 ( Luc vl [19671) . The linear 
limb darkening coefficien ts were interpolated from the tables of 
Wad e & Rucinskil (1 19851) and fixed at Xi(B) = 0.26, x x {V) = 
0.22, xi(I) = 0.165, X] (white) = 0.20. As only monochromatic 
lightcurves are calculated in the lightcurve solution program we 
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Fig. 3. Apparent variation of the atmospheric parameters over 
the phase. Temperature, surface gravity and helium abundance 
were measured from the Gemini/GMOS spectra. The errors are 
statistical errors. 



used the central wavelengths of the band filters. The secondary's 
limb darkening coefficient X2 was varied as the limb darken- 
ing of such heated objects can not be predicted. The tempera- 
ture of the hot component was taken from the spectral analysis 
(r eff ,sdB = 28 400 K). 

For the albedos of the companion values exceeding 1 were nec- 
essary to model the reflection effect. This can be explained with 
processes in the stellar atmosphere that cause a spectral redistri- 
bution of the irradiated energy with wavelength. The third light 




Fig. 4. T e ff - log ^-diagram. The helium main sequence (HeMS) 
and the EHB band (limited by the zero-age EHB, ZAEHB, and 
the terminal-age EHB, TAEHB) are sup erimposed with EHB 
evolutionary tracks from lDorman et all (119931) . The position of 
AS AS 10232 is indicated with a solid diamond. Open diamonds 
mark the position of other HW Vir-like systems (jC harpinet et alj 
2008 : IDrechsel et al.l l200lt iFor et alJl2010t iGeier et al.ll20 1 lb; 
Maxted et al.ll2002t iKlepp & Raucb]l201 lU0stensen et al.lF2 008. 
12010 1: IWood & Safferlll999t lAlmeida et al.ll2012t iBarlow et al.l 
12012b . 



that accounts for the disturbance by a potential third object in 
the system was first varied but then set to as it did not de- 
viate from this value. As the spectral type of the companion is 
very late, radiation pressure does not play a role and the radia- 
tion pressure coefficient fo r the companion was fixed at zero (cf. 
IDrechsel et aB 19951 12001 . 

The rest of the parameters such as the radiation pressure coeffi- 
cient for the primary star, the inclination, the effective tempera- 
ture of the companion and the Roche potentials were adjusted. 
A grid with different fixed mass ratios and different start param- 
eters was calculated. All four lightcurves were analysed simul- 
taneously. As there are so many correlated parameters a unique 
solution can not be found from the lightcurve alone because of 
the degeneracy of mass ratio and radii of the components. In 
Table [3] the parameters of the lightcurve solution with the best 
standard deviation can be found. 

The Euler lightcurves are consistent with this solution as it can 
be seen in Fig. [8] The only difference is a higher albedo of 1.3 
for the companion. As the Euler lightcurves have a good signal- 
to-noise and four different colour bands are covered, the higher 
albedo is probably preferable. 

Fig. [9] shows the equipotential surfaces of both components 
for the case of the best solutions over the phase. It is clearly 
visible that AS AS 10232 is a detached system with a grazing 
eclipse as both components stay inside their Roche lobe. The 
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Fig. 5. Phased SuperWASP lightcurve. The solid line demon- 
strates the best fitting model. In the bottom panel the residuals 
can be seen. The wave pattern that is seen in the residuals can 
be explained by the simplified treatment of the reflection ef- 
fect. Better models of the reflection effect should remove this 
wave pattern. 
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Fig. 6. Phased lightcurves of AS AS 10232 in the B (*), V (x) and 
I (+) band taken from the Carries Hill Observatory: V is shifted 
by 0.05, 1 is shifted by 0.1 



companion is bigger than the primary star and slightly distorted. 
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Fig. 8. B VRI Euler lightcurves with the model from the parame- 
ters of the best solution, but a higher A 2 = 1.3, every lightcurve 
is shifted by 0.05 
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5. System parameters 

As the radial velocity amplitude of the companion is not known 
and a degeneracy in the mass ratio appears in the lightcurve 
solutions we do not have a unique solution of the system yet. 
The solution with the smallest standard deviation from the 
lightcurve model to the measured values is found for a mass ra- 
tio q = 0.34. For this mass ratio we derive an inclination angle 
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o u 



Fig. 9. Equipotential structure of ASAS 10232 corresponding to 
the best fit solution over the phase from 0.5 to 0.5 in steps of 
0.125 
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Table 2. Parameters of ASAS 10232 (TYC 7709-376-1) 



coordinates 


10 23 22 -37 37 00 (J2000.0) 


proper motions 


-28.9 -22.8 [2.4 2.4 35] mas/yr 


mass of the sdB 


M siB 


[Mq] 


0.461 +0.051 


mass of the companion 




[M Q ] 


0.157 ±0.017 


separation 


a 


[*©] 


0.963 ± 0.036 


mean radius of the sdB 


^sdB 


[*0] 


0.179 ±0.011 


mean radius of the comp 


*^comp 


[*©] 


0.256 ± 0.015 


surface gravity (phot) 


logg(sdB) 




5.60 ± 0.02 


surface gravity (spec) 


logg(sdB) 




5.60 ± 0.05 



of i = 65.9° from the lightcurve analysis. Furthermore, the mean 
radii for the components could be determined: R\/a = 0.185 
and Rila = 0.266, where a is the separation of the mass cen- 
tres. Together with the K\ — 81.0 + 3.0kms _1 measured from 
the radial velocity curve and the mass ratio q, the inclination i 
and the period P = 0.13927 d, determined from the lightcurve, 
the masses of the components M\ = 0.461 + 0.051 Mq and 
M2 = 0.157 + 0.017 Mq were derived. Moreover, the sepa- 
ration a = 0.963 + 0.036 R e and therefore the absolute val- 
ues for the radii of the stars R^b = 0.179 + 0.01 1 Rq and 
^comp = 0.256 + 0.015 Rq could be calculated. 
Radius and mass of the companion leads to the conclusion that 
the companion must be an M dwarf. The derived mass for the 
sdB is consistent with theoretical calculation s for the forma- 
tion o f sdBs via the common envelope channel dHan et al . 2002, 
120031) . 

In order to check this solution the surface gravity determined 
from the spectra can be compared to a photometric surface grav- 
ity. This surface gravity is derived via the mass - radius rela- 
tion from the mass and radius calculated from the radial velocity 
curve together with the lightcurve analysis. The surface grav- 
ity determined in this way logg(phot) = 5.60 + 0.02 is in per- 
fect agreement to the spectroscopic surface gravity log g(spec) = 
5.60 + 0.05. 

6. Conclusions 

We discovered the new eclipsing sdB star ASAS 102322-3737.0, 
which shows a grazing eclipse and a huge reflection effect that 
is the only contribution of the companion to the optical light of 
the system. It has an orbital period of P = 0.13927 d and an 
inclination of i = 65.9°. A unique solution can not be found as 
there exists a degeneracy in the mass ratio. The best lightcurve 
solution was calculated for a mass ratio of q — 0.34. The masses 
of the components for this solution M s[ jb = 0.461 + 0.051 Mq 
and M comp = 0.157 + 0.017 Mq are typical for an sdB and a 
late main sequence star. Also the spectroscopic and the photo- 
metric surface gravity are in agreement. As ASAS 10232 is very 
bright it is possible to apply rapid high resolution time series 
spectroscopy to search for spectral features and the RV semi- 
amplitude from the companion to resolve the degener acy in the 
mass ratio, as done for AA Dor (Vuckovic et al. 2008). 
Some of the sdBs were found to be pulsating. There are two 
classes of sdB pulsators, p-mode pulsators with periods of some 
minutes and g-mode pulsators with periods of 0.5 to 2 hours, 
which are separated by their atmospheric parameters. Until now 
only one pulsating HW Virginis system was found, NY Virginis 
dVuckovic et alj|2007l) . Such systems are very interesting as it 
is possible to compare the results from asteroseismology to the 
results from the lightcurve analysis of the eclipses. If it is pulsat- 
ing, the parameters of ASAS 10232 would suggest it to be most 
likely a g-mode pulsator. However, the Euler lightcurve with the 



best signal-to-noise shows no sign of pulsation. But the strong 
reflection effect could hide the pulsations. 
Period variations have been found for almost all of the HW Vir 
systems that have accurate eclipse timings covering more than 
five years, which may be attributed to the presence of a third 
body. In several cases the third bodies are likely to be one (or 
two) giant planet(s). These discoveries were unexpected, be- 
cause it is regarded difficult for giant planets to form around 
main sequence binaries due to the short life time of circumbi- 
nary disks. In addition, such planets may not be able to sur- 
vive common envelope (CE) evolution. Instead, it has been sug- 
gested that these circumbinar y planets are of second generation 
(IZorotovic & S chreiber 2012) , that is formed from the instability 
of a post-CE disk. |Zorotovic & Schreiberl (1201 2b proposed also 
an alternative scenario for the period variations due to processes 
acting in deeply convective secondary stars. 

HW Vir has b een monitored for more than 28 years now and 
Lee et al .] (120091) found two sinusoidal variations of the light- 
travel time for HW Vir from 24 years of data, suggesting the 
presence of two substellar objects orbiti ng the close binary. New 
observations by iBeuermann et al.| (120121). howeve r, deviate sig- 
nificantly from the prediction of iLee et al.1 (120091) . The new so- 
lution involves one planet an d a brown dwarf or low -mass star 
orbiting the HW Vir binary (IBeuermann et al.ll2012l) . The new 
solutio n was found to be stable in contrast to that of ILee et"aT1 
(120091 

These findings suggest that the probablity to find period vari- 
ations in ASAS 10232 as well is high, whether due to substellar 
objects or the active secondary. The lesson learned from HW 
Vir is that long term and dense monitoring is a prerequisite 
(Beuermann et al. 2012). ASAS is the third brightest HW Vir 
system known, only one magnitude fainter than HW Vir, which 
allows us to use readily accessible small telescopes to obtain 
lightcurves. Hence it is a promising target for amateur and high- 
school observatories t o team up with professional astronomers 
dBackhaus et al]l2012l) . 
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